Optical A transitions and quantum computing in the 15 N-V Center in Diamond 
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We present a thorough derivation of the excited state energy levels of the negatively charged 15 N-V~ center 
in diamond for the strong applied electric field case. We show that in the 15 N-V~ center a spin non-conserving 
two-photon A transition exists that is mediated by the hyperfine interaction, which provides the possibility to 
write quantum information. Using second order perturbation theory we obtain a A transition rate of the order 
of 10 MHz at room temperature, which allows for approximately 10 4 quantum logic operations within the spin 
coherence time r d (T = 300 K) « 1ms of the 15 N-V" center. 
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Recently tremendous attention has been drawn by single 
Nitrogen-vacancy defect centers (N-V centers) in diamond 
due to their promising properties for the experimental real- 
ization of quantum computation (TJ. In particular, the N-V 
center in diamond is well suited for studying electronic and 
nuclear spin phenomena, since its spin can be both initialized 
and read out via a stable optical spin conserving transitions 
10 [3]. The N-V center appears in two forms, neutral and neg- 
atively charged, and it is the negatively charged center (N- 
V~) that we are interested with in this rapid communication. 
The N-V~ stands out among solid-state systems because its 
electronic spin can be efficiently prepared, manipulated and 
measured with optical and microwave excitations. Recent ex- 
periments have conclusively demonstrated that the N-V center 
in ultrapure diamond lattice shows the longest room tempera- 
ture spin dephasing time ever observed in solid state systems 
(T2=1.8 ms) J4). The electron spin coherence lifetime in 
the N-V~ center is limited by its hyperfine interaction with 
the Carbon lattice in the diamond structure [5]. At very low 
concentration of paramagnetic impurities the loss of electron 
spin coherence is minimized and spin quantum effects could 
be observed at room temperature over hundred of nanome- 
ters, these results open the door to coherent manipulation of 
individual isolated nuclear spins in a solid state environment 
at room temperature with potential applications to quantum 
computation, therefore the understanding of the hyperfine in- 
teraction in this solid state system is of great interest. 
In this rapid communication we show how the hyperfine inter- 
action gives rise to a spin non-conserving two-photon A tran- 
sition with a rate of the order of 10 MHz at room temperature, 
which allows for approximately 10 4 optically controlled quan- 
tum logic operations within the spin coherence time at room 
temperature. The A transition can be understood by consider- 
ing the hyperfine interaction of the electron's spin in the N-V~ 
center with the spin- 1/2 nucleus of the isotope Nitrogen 15 N. 
A quantitative result for the number of quantum logic oper- 
ations is obtained by using second order perturbation theory 
and numerical diagonalization of the Hamiltonian describing 
the excited energy levels of the 15 N-V~ center when an ex- 
ternal electric field is applied. Our results clarify the origin of 
the two-photon transitions that have been experimentally ob- 



served in Refs. 0121. 

The N-V center in diamond consists of a substitutional nitro- 
gen atom next to a Carbon vacancy giving a centre with C3„ 
symmetry. In the negatively charged state, the extra electron 
is located at the vacancy site forming a spin 5=1 pair with 
one of the vacancy electrons. The energy level structure of the 
N-V~ center was established by combining optical electron 
paramagnetic resonance and theoretical results. Lenef and 
Rand earlier work gave a group theoretical description of the 
molecular orbitals and state configuration of the N-V~ center 
by using linear combinations of sp 3 orbitals [8|. Recent nu- 
merical ab intio calculations conducted by Hossain et al. in- 
dicated the following orthonormal orbitals which irreducibly 
transform according to the C3,; symmetry group : u = Sjy, 
v = (si+s 2 +s 3 -3f3s N )/V3Si, e x = (2s 3 -si-s 2 )/V3S 2 , 
e y = (si — s 2 )/S 2 , where sn, si, s 2 and S3 denote sp 3 or- 
bitals pointing to the Nitrogen and three Carbons neighbor- 
ing the vacancy, Si = ^l + 2a-3/3 2 , S 2 = \/2 - 2a, 
a = (si\s 2 ) and (3 — (s 3 |u) are normalization constants [9|. 
The energy levels are traditionally labeled according to the 
orbital and spin multiplicity, which is related to the number 
of states having the same energy. For the negatively charged 
N-V center, the ground state is a spin triplet state 3 A2, with 
a zero-field splitting of D gs — 2.87 GHz between spin sub- 
levels m s = and m s = ±1 [see Fig. [JJ. 
The excited state 3 E is also a spin triplet, associated with 
a broadband photoluminescence emission with zero phonon 
line at 1.945 eV, which allows optical detection of single N-V 
defects using confocal microscopy. It is now well established 
that a metastable state ^1 is place between the ground and 
excited triplet states [10]. Although an extensive documenta- 
tion have been made to unveil the excited state structure of the 
N-V center from first principles, it remains a topic of current 
research due to its possible quantum applications [11]. 
The Hamiltonian suitable to describe the 3 E excited state of 
the negatively charged 15 N-V~ center is given by lfT2l 

H=H el +\L-S+p[(L-S) 2 +l-L-S-l-L(L+l)S(S+l)}+AI-S. 

(1) 

The first term in Eq. ([TJ is the dominant term in the Hamil- 
tonian and is given by H e i = F ■ D, where F is the external 
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FIG. 1: Schematic diagram for the energy levels of an NV~ center 
depicting the spin triplets for the ground state 3 A2 and the excited 
state 3 E. The splitting due to spin orbit interaction is shown for the 
excited state. 



applied electric field in the N-V defect and D is the dipole 
moment operator, respectively. From earlier work done by 
Hossain et al. we know that the only nonzero dipole mo- 
ments components are D Ca . and D Ey with equal dipole mag- 
nitudes, i.e. |D e J = |D e | = D, respectively. The dipole 
moment transverse terms cause a linear splitting of the degen- 
erate 3 E state levels given by ±DyJ(F x + f x ) 2 + (F y + f y ) 2 
lfT2l . where we have taken into account that the presence of 
strain has the same effect as a weak electric field denoted by f . 
The second and third terms in Eq. ([T]i describe the spin-orbit 
and the spin-spin interaction in the LS coupling approxima- 
tion fl3l . respectively. 

With this information and restricting ourselves to the 
mi = ±1 state sublevels of the spin triplet we can write 
the electric field, the spin-orbit and spin-spin interactions in 
the following matrix form 
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where we have chosen the basis \mim s ) in the following or- 
der [| - 11), | - 1 - 1), | - 10), 1 10), |11), 1 1 — 1>] for numerical 
purposes and we have assumed that the electric field is applied 
in the y direction. The direct product between matrices takes 
into account the nuclear spin subspace, i.e. mi — ±1 /2. 
One can see that the spin-orbit interaction does not couple 
terms with different orbital projections and produces just an 
energy splitting between the states m s = ±1. The spin- 
spin interaction however couples the first and the last state 
with different orbital projections and produces an energy split- 
ting between the states m/ = ±1. The Hamiltonian given 
in Eq. Q can be solved exactly and it is easy to convince 
oneself from the exact solution that there is no overlapping 
between eigenstates with different spin projections. Hence, 
the spin-orbit and spin-spin interaction cannot account for the 
mixing of the excited states and therefore cannot give rise 
to spin non-conserving transitions. If we want to allow for 
spin non-conserving transitions between one ground state and 
another via a common excited state we must introduce the 
nuclear-spin interaction. The hyperfine interaction will mix 
the eigenstates with different spin projections allowing us to 
have a transition to two different ground states. This particu- 
lar transition is a second order process that can be calculated 
by using Fermi's golden rule and is commonly known as a A 
transition. Therefore, the last term in the total Hamiltonian 
given in Eq. ([T]) is the nuclear-spin interaction and this term 
will be responsible for the mixing of the excited states, which 



will give rise to A transitions in the 15 N-V~ center in dia- 
mond. The hyperfine interaction can be written in the above 
basis as 



H ns =A 
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(3) 



where A is the excited state hyperfine coupling constant which 
for the 15 N nuclear spin case is A m 60 MHz [14|. The 
eigenvalues of the total Hamiltonian given by Eq. ([T]) are 
plotted as a function of the electric field in Fig. [2] The en- 
ergy diagram in Figure|2]shows the mixing between the states 
rn s = ±1. This mixing will allow us to induce A transitions 
between these states in order to write quantum information by 
using circularly polarized light. The 15 N-V~ center has opti- 
cal transitions between the ground state 3 A2 and the excited 
state 3 E at zero magnetic field that are spin conserving and 
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FIG. 2: Energy diagram showing the eigenvalues of the total Hamil- 
tonian given in Eq. {I} as a function of the applied external elec- 
tric field. The lower branch shows the mixing between the states 
m a = ±1. The eigenvalues are doubly degenerate. 



will allow the readout of information through the electron's 
spin lfl2ll . In order to write information we need to have op- 
tical transitions from two ground state levels via a common 
excited state. The induced transition probabilities for the case 
of circularly polarized light around the z-axis is denoted by 
<t and is proportional to the modulus squared of the matrix 
element given by (i/'.j |<r ± where ipi and ipj are electronic 
wave functions for the initial and final state, respectively, and 
(7 = ^^Ji„ ± iy n )/y/2, where the summation are over 
the electrons in the N-V~ center. 

Suppose that in the absence of light the 15 N-V~ system is 
prepared in the ground state \mim s )\mi) = |00)| |). If a 
strong external electric field with magnitude F is applied and 
the crossing between different states occur at a given energy 
e„ for the n th excited state as illustrated in Fig. [5] then the 
ground state will be coupled to an excited state given by 



E 
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The coherent superposition of states given in Eq. Q allows 
the possibility to have optical transitions between two ground 
states with different spin projections, i.e. spin non-conserving 
transitions, by using either a right or left circularly polarized 
light. In order for the final total angular momentum to equal 
the initial angular momentum, the nuclar spin projection flips 
[See Fig. [3j. The amplitude A± for the A transition in the 
dipole approximation from the initial state 1 00) | |) to the final 
state |01)| i) through either right or left circularly polarized 
light is given in second order perturbation theory by 
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where e is the electron charge, £ is the amplitude of the drive 
optical field, eo is the energy of the initial ground state and 



u a ± is the frequency of the circularly right or left polarized 
light, respectively. Writing the explicit form of \^„) in terms 
of the basis vectors in Eq. |5]l we get 
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The value of the CJJ, imsm; coeffi- 



cients and eigenvalues e n can be obtained by numerical diag- 
onalization of the total Hamiltonian for a given value of the 
applied external electric field Fq, [see Table [[). In the subse- 
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0.5633 


0.0056 


-i0.0021 
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-0.1011 
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-i0.0004 
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0.5754 
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0.1344 
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12 


10.9909 


-0.0271 


-0.6961 


-i0.6962 


i0.0122 



TABLE I: Numerical values for the coefficients C^ l[ msrni and eigen- 
values e„ of the excited state levels that contribute to the transition 
amplitude A± for the case F = 1.85 MV/m fl5l . 

quents calculations we will work with right circularly polar- 
ized only. In order to calculate the transition rate between the 
initial ground state and the final ground state we use Fermi's 
golden rule 
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where p(uj) represents the allowed photon frequencies and can 
be written in its simplest form as a Lorentzian function 
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where e is the energy of the final ground state and 7 is the 
decoherence rate. For the resonant case, i.e. ui a + — Lo' a+ — 

(e' - e ), Fermi's golden rule is given by W^|oo)|T)^|oi)|i) = 
4\A\ 2 /'jh 4 . Taking the value of the dipole moment matrix 
elements as eo^-i = ecr^ = \d\ in Eq. (|6j) we get 
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where A„ = ui a + — (e„ — eo) is the detuning and we have used 
the relation £\d\ = Ml r where is the Rabi frequency. 
In order for perturbation theory to be valid we require that 
|A„| ^ i7/j|C™ *| |, hence we have to tune the frequency Lo a + 
to achieve this. Assuming that the detuning is taken in such 
a way that |A„| ps 30rii?|C" *| |, where the Rabi frequency is 
a few GHz [16|, then the amplitude is \A\ 2 « lO- 27 h 4 tt 4 R . 
Using the spin coherence time of j^ 1 w 1ms flU, we can 
estimate the A transition rate to be of the order of 

W|oo>|T>->|oi}|j.) ~ 10 MHz. (10) 

Based on the obtained transition rate and for coherence life- 
times of around the milliseconds [4|, we estimate that the 15 N- 
V~ center in diamond should be capable of around 10 4 quan- 
tum operations per coherence lifetime at room temperature, 
which confirms that this solid state system is a very good can- 
didate for quantum computing flTl . 

Suppose now that we can prepare the 15 N-V~ center in such 
a way that all the spin population is in the ground state level 
rn s = |0). We have shown that an optically driven A transition 
will change the nuclear spin projection because of the spin 
non-conserving transition. Therefore if we optically pump 
the ground state level so that we have a A transition of the 
form \m s = 0) — > \m s = ±1), then the ground state pop- 
ulation will transfer according to the two possible A transi- 
tions \m s = 0>| t>-> \m a = ±1>| |> or \m s = 0)| 
\m s — ±1)| f>, with flipped nuclear spin projection [see 
Fig. [3). Note that both A transitions are completely indepen- 
dent from each other. Then, by means of resonance condition, 
it is possible to produce an arbitrary superposition of the form 
\^w) = ot\m s = 0) + f3\m s = 1), which represents the qubit. 




m t = 1 0> m s = 1 0> m, = | J, > m 1 =|o>m s =|o>m,=|t>,e 



FIG. 3: The A transition optically driven with (j + -polarized laser 
beams will flip the nuclear spin projection in order for total angular 
momentum to be conserved. The A transition will equally popu- 
late both ground state levels due to the two-fold degeneracy between 
them. 



In order to read this superposition state the electric field 
F y needs to be switched off, leading to A = 0. This means 
that for F y — and weak strain the two-photon transitions 
are spin conserving, which allows for the detection e.g. 



| a | 2 through the linear response to a <r + -polarized laser 
beam. Single-qubit operations of the qubit state \^f w ) can 
be performed by means of the spin conserving and spin 
non-conserving transitions, which give rise to phase and 
amplitude shifts, respectively. 

In conclusion, we have demonstrated that there exists a A 
transition in 15 N-V~ center in diamond that is mediated by 
the hyperfine interaction which allows us to write quantum 
information in this solid state system. Readout can be 
performed when the electric field is tuned to spin conserving 
transitions. Using second order perturbation theory we 
have obtained a A transition rate of the order of 10 MHz 
at room temperature, which allows for approximately 10 4 
quantum logic operations within the spin coherence time 
r d (T = 300 K) « 1ms of the 15 N-V~ center. Our findings 
clarify the origin of the two-photon transitions observed 
in Refs. flU [71, which is essential for the development of 
a scalable quantum network made of 15 N-V~ centers in 
diamond, given the possibility to produce a photonic crystal 
structure in diamond |fT8l[T9l . 
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